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TECHNICAL PAPER 

A NEW STATE RECONSTRUCTOR FOR DIGITAL CONTROL SYSTEMS 
USING WEIGHTED-AVERAGED MEASUREMENTS 

I. INTRODUCTION 


Consider the usual linear continuous-time plant driven by a zero-order-hold with sampled output, 
as shown in Figure 1 . The most common approach to reconstructing the state of this system is undoubted- 
ly the state observer [1], However, it is well known that the state observer has some undesirable 
characteristics. Firstly, it is a dynamical system in itself, and, hence, adds additional states and eigen- 
values to the system. Secondly, as a consequence, the reconstructed state is normally an approximation 
to the true state and is usually not a good one early in the state reconstruction process, unless the initial 
state of the system is well known. 


U(kT), 



Figure 1 . Linear continuous-time plant driven by a zero-order-hold with a sampled output. 

The work in References 2 to 7 focuses on a relatively new state reconstructor that has neither of 
these problems. If the plant is noise-free and can be modeled precisely, the output of this state reconstruc- 
tor exactly equals the true state of the plant and accomplishes this without any knowledge of the plant’s 
initial state. Besides, it adds no new states or eigenvalues to the plant equation in the discrete time model 
of the system. In fact, it does not affect the plant equation in the discrete time model at all; it affects the 
measurement equation only. It is characterized by the fact that discrete measurements are generated every 
T/N seconds and input into a multi-input/multi-output moving average process. The output of this 
process is sampled every T seconds and utilized in reconstructing the state of the system. Should process 
noise, measurement noise, or modeling errors be a problem, it can be catenated with an observer or a 
Kalman filter for a synergistic effect. 

In this paper, another new state reconstructor is presented. This one has the same advantages, 
i compared to the state observer, as the one developed in References 2 to 7. However, it is unique in the 

i way it achieves state reconstruction. A continuous-time output vector from the plant is convoluted with a 

[ weighting-function matrix whose elements are time dependent. This result is integrated over T second 
intervals to generate weighted-averaged measurements, every T seconds, that are used in the state recon- 
struction process. In actuality, the weighting-function matrix is the impulse-response matrix for a multi- 
input/multi-output system and the weighted-averaged measurements are the forced responses of this 

> 

i 

1 


/ 




system to the plant output over T second intervals. As a result, an equally valid way of generating the 
weighted-averaged measurements is as follows. The continuous-time output vector from the plant is 
input into a continuous-time multi-input/multi-output system that is initially at rest. After T seconds, the 
output of this system is the weighted-averaged measurement vector. This process is repeated every T 
seconds. 

This new state reconstructor is presented in the remainder of this paper. In Section II, some pre- 
liminary mathematical formulations pertinent to the development of it are derived. In Section III, the new 
state reconstructor is presented. An example of it is given in Section IV. Some final comments about it 
and its implementation are made in Section V. 


II. PRELIMINARY 


For the plant in Figure 1, 2i(t)€R n is the state, jj(kT)€R r is the control input, y s (kT)€R m is the 
standard output or measurement vector, F6R nxn is the system matrix, G€R nxr is the control matrix, and 
C s €R mxn is the standard output matrix. It is well known that this system can be modeled at the sampling 
instants kT by the discrete state equations [8] 


_x(k + 1 ) = A_x(k) + Bju(k) 

(1) 

y s (k) = C s x_(k) , 

(2) 


where k is the usual shorthand notation for time kT, 


4>(t) = [(sI-F) 1 ] , (3) 

A = <|>(T) , (4) 


and 



<J>(t)€R nxn is the state transition matrix. A€R nxn is the system matrix and B€R nxr is the control matrix for 
the discrete state equations (1) and (2). 
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Now consider the plant in Figure 2, which is a generalization of the one in Figure 1 . In addition to 
the standard output y s (kT), it has the output y^ A (kT)eR 4 found by differencing the weighted-averaged 
measurement y A (kT)€R q with the delayed-input term Hjj[(k-1)T] where H€R qxr . The vector y^kT) is 
generated by convoluting the continuous-time output z(t)€R p for (k-l)T < t kT with the weighting- 
function matrix 0(t)€R qxp and integrating this result over T second intervals. Observe that this is equiv- 
alent to inputting ^(t) for (k — 1 )T < t =£ kT into the multi-input/multi-output system defined by the trans- 
fer-function matrix 


0(S) = f 0(t) e st dt € R qxp 

o 

where s is the Laplace transform variable. If the system is at rest at time (k— 1 )T, its output at time kT will 
be y A (kT). It follows from the definitions of ju(t) and z.(t) that the output matrix C A €R pxn . 

For the system in Figure 2, equations 1 to 5 apply. Observe that 


y r (kT) 


’ys(kT) " 
_y A '(kT)_ 


where 


( 6 ) 


y A '(kT) = y A (kT) - Hu[(k-1)T] 


(7) 


and 


.kT 


kT 


y A (kT) = f 0(kT-€) z(0 d£ = f 0(kT— £) C A x(£) d£ . 
(k-l)T (k-l)T 


( 8 ) 


It is known that [8] 

x(£) = 4»(^— 1 0 ) x.(t 0 ) + / 4 >(£-t) Gu(t) dT . (9) 

to 

Letting t„ = kT in equation (9) and substituting this result into equation (8) yields 

• 

( 10 ) 


kT I 

r kT 


J 0(kT-£) C A <f>(e-kT) d£ 

x(kT) + J 0(kT-£j) C A 

J t) Gu(t) dT 

(k-l)T 

(k-l)T 

kT 
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Since _u(t) = u[(k-l)T] for (k-l)T =£ t < kT, equation (10) can be written as 


y A (kT) = Dx.(kT) + Eu[(k-1)T] 


(H) 


where 


D = f 0(kT-O C A cJ>(e-kT) d£ 

(k-l)T 


( 12 ) 


kT 


E = f 6(kT-0C, 

(k— I )T 


/ <t>(^— t) dT 


kT 


G 


the transformation X = £-kT in equation (12) produces 


r 

D = J 6(-X) C A A(\) dX 


(13) 


(14) 


A(X) = (t>(X) 


(15) 


he transformations X = £-t and p .= £-kT successively in equation (13) yields 

r° 

E = J 6(-p) C A B(p) dp 

-T 
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From equations (2), (6), (7), and (11), 


y T (kT) = 


Cs 


' 0 “ 


_x(kT) + 

E-H 

D 



u[(k-l)T] 


Let 

f ° 

H = E = J 0(-p) C A B(p) dp . 

-T 

In this case, the discrete state equations for the system in Figure 2 become 
_x(k+l) - A_x(k) + Bufk) 
y T (k) = C T _x(k) 

where the output matrix C T €R ,m+t|,xn is 


Ct — 


Cs 

LD 


This follows from equations (1), (16), and (17). 

III. THE NEW STATE RECONSTRUCTOR 


(16) 


(17) 


(IB) 

(19) 


( 20 ) 


A block diagram of the new state reconstructor catenated with the plant is shown in Figure 3. Its 
output is y-r'(kT)€R n . Let the matrix Cx€R nx<m + q) be the pseudo-inverse of Ct- That is, let 

Ct = (C t t CtF 1 C t t , 

as shown in Figure 3, where Cj is defined by equation (20) and (Ct t Cj) is assumed nonsingular. Then 
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( 21 ) 


y-r'(kT) = C*yj(kT) = (C T T C T Y l C T T y T (kT) 


It was shown in Section II that the discrete state equations for the system in Figure 3 from u(kT) to yj(kT) 
are given by equations (18) and (19). Using these and equation (21), the discrete state equations for the 
system in Figure 3 become 


x(k + 1 ) = A x(k) + B u(k) 


y T '(k) = x(k) 


( 22 ) 


Hence, the state reconstructor exactly reconstructs the state of the system and does so without affecting 
the plant equation in the discrete time model of the system. 

It was assumed that (C T T C T ) is nonsingular. This places certain restrictions on C T €R <m + q)xn [9]. 
First, it is necessary that (m + q) 2= n, or equivalently, q 2 s (n-m). This is readily satisfied since q, the 
number of rows in 0(t)€R qxp , can be chosen arbitrarily. Secondly, it is necessary that rank (C T ) = n. By 
virtue of equations (14), (15), and (20), rank (C T ) is affected by <j>(t), 0(t), C s , and C A . In most problems: 
4>(t) will be fixed; 0(t) can be chosen arbitrarily; C s and C A will lie somewhere in between. Hence, given 
4>(t), if 0(t) can be chosen, and perhaps C s and C A to some extent, so that rank (C T ) = n, the state of the 
plant can be reconstructed according to equation (22) for k s* 1 . 


IV. AN EXAMPLE 


Consider the double integrator plant driven by a zero-order-hold as shown in Figure 4. The output 
z(t) = x | (t) is to be used in generating the weighted-averaged measurements for state reconstruction; 
there are no standard measurements in this example. Manipulating this plant into the format of Figure 3 
yields 


F = 


0 1 
0 0 


G = 


O' 

1 


C A = [10], 
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Figure 4. Double integrator plant in the example. 


where Cs is a null matrix. It is straightforward to show that 


A(t) = 4>(t) = 


"1 t" 
_0 1 _ 


f 4>(\) dX = 2 

0 lo t . 


and 



p _ 


r-t 2 i 

B(t) = 

/'(f)(X) dX 

G = 

2 1 


0 


. t . 


It is clear that n = 2, m = 0, and p = 1 . The state reconstructor in Figure 3 requires that q 2^ (n-m) or 
q 2. Choose q = 2. Hence, the weighting-function matrix 0(t)€R qxp becomes 0(t)€R 2xl . Choose 



"e,(t>“ 


1 

0(t) = 


= 



_ 0 2 (t) _ 


L e -(2/T)t 

-T 


Taking the Laplace transform of 0(t) produces 


0(S) 


re.(S) 


L0 2 (S)J 



(23) 
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Observe that 0|(S) is an integrator and 02(S) is a low-pass filter with a time constant of T/2 seconds. 
Hence, in implementing the state reconstructor for this example, the weighted-averaged measurements 
could be generated in either of two equally valid ways. One way is to determine the convolution integral 
shown in Figure 3 using equation (23). The other is to input z(t), (k-l)T<t^kT, into the system 0(S) 
with zero initial conditions at t = (k-l)T. At t = kT, the output of this system is y A (kT). See how 0(S) 
can affect measurement noise in z(t) if it is present. Choosing 0(t), in general, to reduce the effects of 
measurement noise is an area for future study. 


Having chosen 0(t) in this example, the other parameters in the state reconstructor can be found 
using the formulations in Section III. Doing this, it turns out that 


|^(1 - e 2 ) 5 (3e " 2 - 1} 

'($ & 
m ( 3 . ' 


r i! 

6 


The new state reconstructor is now completely defined for this example. 


V. FINAL COMMENTS 


A new state reconstructor for a linear continuous-time plant driven by a zero-order-hold has been 
presented. It is unique in that it generates weighted-averaged measurements from a continuous-time 
output and utilizes these in reconstructing the state of the plant every T seconds. It has the same advan- 
tages, compared to the state observer, as the state reconstructor developed in References 2 to 7. These are 
enumerated in Section I. Should process noise, measurement noise, or modeling errors pose a problem 
for it, it could be catenated with an observer or a Kalman filter for a synergistic effect. Determining the 
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effects of noise and modeling errors on the new state reconstructor is an area for future study. Reducing 
their effects by judiciously choosing the time-dependent elements of the weighting-function matrix is 
another area for future study. 

With regard to implementing this new state reconstructor, it seems natural to generate a vector of 
weighted-averaged measurements wherever a scalar continuous-time output exists. Every T seconds, all 
measurements, weighted-averaged and standard, could be sent to a central processor. There, the 
weighted-averaged measurements could be properly combined and all further calculations and operations 
performed to reconstruct the state of the system. There are several ways of calculating the weighted- 
averaged measurements. They could be done in microprocessors using digital approximations or else in 
analog circuitry. In either case, a convolution integral, transfer function, or state equation approach is 
possible. This is also an area for future study. 
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